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Specific features of the equilibrium current-carrying state of a Josephson tunnel junction between 
diffusive superconductors (with the electron mean free path / smaller than the coherence length fo) 
are studied theoretically in the ID geometry when the current does not spread in the junction banks. 
It is shown that the concept of weak link with the jump $ ~ 1 of the order parameter phase exists 
only for a low transmissivity of the barrier F < Z/£o <C 1- Otherwise, the presence of the tunnel 
junction virtually does not affect the distributions of the order parameter modulus and phase. It 
is found that the Josephson current induces localized states of electron excitations in the vicinity 
of the tunnel barrier, which are a continuous analog of Andreev levels in a ballistic junction. The 
depth of the corresponding "potential well" is much greater than the separation between an Andreev 
level and the continuous energy spectrum boundary for the same transmissivity of the barrier. In 
contrast to a ballistic junction in which the Josephson current is transported completely by localized 
excitations, the contribution to current in a diffusive junction comes from whole spectral region near 
the energy gap boundary, where the density of states differs considerably from its unperturbed value. 
The correction to the Josephson current j'($) in the second order of the barrier transmissivity, 
which contains the second harmonic of the phase jump $, is calculated and it is found that the 
true expansion parameter of the perturbation theory for a diffusive junction is not the tunneling 
probability F, but a much larger parameter W = (3£o/4Z)r. This simplifies the conditions for the 
experimental observation of higher harmonics of j(<&) in junctions with controllable transmissivity 
of the barrier. 

Fiz. Nizk. Temp. 25, 230-239 (March 1999) 



X 



I. INTRODUCTION 

In recent years, considerable advances have been made 
in technology of preparing low-resistance tunnel junc- 
tions with a comparatively high barrier transmissivity 
(tunneling probability) T. This primarily applies to con- 
trolled break-junctions Q| as well as systems based on 2D 
electron gas Q , whose conductivity undergoes a crossover 
from tunnel to metal type upon a change in the barrier 
parameters. The problem of calculation of the Josephson 
current through a junction with an arbitrary transmis- 
sivity in the ballistic regime (with the electron mean free 
path I much greater than the coherence length £o) was 
solved by many authors || on the basis of the model of a 
single-mode junction with current-carrying banks ensur- 
ing a rapid "spreading" of supercurrent and the equality 
of the order parameter modulus A near the barrier to its 
bulk value (the "rigidity" condition for A). 

In the ID geometry (e.g., a planar junction or a super- 
conducting channel with a tunnel barrier Q), the prob- 
lem is complicated considerably due to the change in the 
order parameter and the quasiparticle energy spectrum in 
the vicinity of the junction, which makes a contribution 
to the phase dependence of the current Antsygina 
and Svidzinskii Q determined the corresponding correc- 
tions to j($) of the order of T 2 for a pure (I S> £o) su- 
perconductor in the limit of low transmissivity f< 1: 



5j($) = -a(T)7(A)r f sin - sin 2$ 



a(T)~l, (1) 



1(A) = (ir/4)ev F v F TA = 4(A)/ tanh(A/2T), (2) 

where vp is the density of states, vf the Fermi velocity, 
and I C {A) the critical current through the junction. 

In a diffusive superconductor (the "dirty" limit I <C 
£o = \/D/2A, D = vfI/3 is the diffusion constant), the 
calculation of the Josephson current for an arbitrary T is 
hardly possible [|| even in a simple model disregarding 
the variation of the order parameter in the vicinity of the 
junction. As a matter of fact, the boundary conditions 
for isotropic Green's functions g(r,tit 2 ) at the junction, 
obtained by Kupriyanov and Lukichev S , 



l{gVg) L = -l(mU = \ (^) [h,9R] 



(3) 



where d(cos 9) is the tunneling probability for an elec- 
tron impinging the barrier at an angle 9, and the sub- 
scripts R and L mark the value to the right and left of 
the barrier, are valid only within the first order in small 
angle-averaged transmissivity T — (fj,d(fi)). Lambert et 
al. |n| proved that the derivation of the boundary con- 
ditions in the general case (d < 1) is reduced to an anal- 
ysis of a system of nonlinear integral equations for the 
terms in the expansion of the averaged Green's function 
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9{ r iP) — 9( r ) +P9i{ r ) + ■ ■ ■ over Legendre polynomials. 
This problem can be solved only for T <C 1 by expanding 
the right-hand side of Eq. (3) into a power series in T, 
which was used in Ref. for calculating the corrections 
to the Josephson current of the order of T 2 . 

In this paper, we pay attention, first of all, to the fact 
that the problem of calculation of the current-phase re- 
lation for a diffusive junction in the ID geometry has the 
sense only in the case of low transmissivity of the bar- 
rier. Indeed, simple estimates obtained on the basis of 
the well-known formula for j(G>) in the first order in T, 

j ($) = /(A)tanh(A/2T)sin$ (4) 

(which coincides, according to the Anderson theorem, 
with the Ambegaokar-Baratoff result for a pure super- 
conductor jll) ) , show that even for small V ~ l/£ <^ 1 
the critical current through the junction becomes of the 
order of the bulk thermodynamic critical current n s ev sc , 
where v sc ~ l/m£o is the critical velocity of the con- 
densate, n s ~ mi/pDA its density, m the electron mass 
(h = 1). Thus, for r 3> Z/£o the tunnel junction does 
not play any longer the role of "weak link" with the 
jump of the order parameter phase $ and other fea- 
tures of a Josephson element. This follows even from 
the boundary conditions Eq. (3) if we use the estimate 
~ g/£o in the vicinity of the junction, which leads 
to \g L ,9R] ~ sin$ ~ £o/ff « 1 for T » Z/£ @- This 
criterion of weak link can be also formulated in terms of 
the conductance of the system in the normal state: the 
resistance of the junction must exceed the resistance of a 
metal layer of thickness £o- 

From this it follows that the parameter 

W = (3£ /4J)r » r (5) 

plays a fundamental role in the theory of Josephson effect 
for diffusive junction (the factor 3/4 is chosen for conve- 
nience of notation) . We can attach to this parameter the 
meaning of the effective tunneling probability for Cooper 
pairs, which is higher than the conventional probability 
r of quasiparticle tunneling. Small values of W <C 1 cor- 
respond to "weak link" conditions (Josephson effect); for 
W > 1, the presence of a tunnel barrier virtually does not 
affect the supercurrent flow and the distribution of the 
order parameter in a diffusive superconductor. Moreover, 
we can expect that just W and not T is a true parame- 
ter of the expansion of j(<&) in the barrier transmissivity. 
Indeed, the dependence of the Josephson current on the 
mean free path is absent only within the main approx- 
imation in r, Eq. (4) and, therefore, it must be mani- 
fested in higher-order terms of the expansion of j(3>) in 
the emergence of additional dimensionless parameter £ /Z 
in them, which vanishes at I — > oo. An analysis of correc- 
tions to the current-phase dependence of Eq. (4) , carried 
out in Sec. 4 of this article in the next order in W, con- 
firms these considerations and proves that the corrections 



~ r 2 to the Josephson current obtained in Ref. and 
associated with the corrections to boundary conditions 
Eq. (3), are much smaller and insignificant in fact. 

Another important result of the analysis of the current- 
carrying state of a diffusive Josephson junction is the 
conclusion concerning the emergence of localized states 
of electron excitations in the vicinity of the barrier. This 
phenomenon is well known for a ballistic tunnel junction 
fll~3| , [i4j in which discrete energy levels 

e„($) = ±A(1 -dsin 2 $/2) 1 / 2 , (6) 

associated with Andreev localization of electron excita- 
tions near the jump in the order parameter phase, split 
from the continuous spectrum in the current-carrying 
state. A similar phenomenon also takes place in a diffu- 
sive junction in which, however, isolated coherent energy 
levels cannot exist due to electron scattering at impuri- 
ties and defects. In this case, the most adequate descrip- 
tion of the variation of the energy spectrum of excita- 
tions is the deformation of their local density of states 
N(e,r) = Re u R (e,r) (u R is the diagonal component 
of the retarded Green's function for the superconductor) 
which is assumed for brevity to be normalized to its value 
vp in the normal metal. In the absence of current, the 
density of states in a homogeneous superconductor has 
the standard form N (e) = |e|9(e 2 -A 2 )/Ve 2 - A 2 (6(x) 
is the Heaviside function) with root singularities at the 
gap boundaries. In the current state, the momentum p s 
of the superfluid condensate plays the role of a depairing 
factor smoothing the singularities of N(e) and reducing 
the energy gap 2e* by A - e*(p s ) ~ (Dp 2 ) 2 / 3 Q. In 
the vicinity of a weak link, a similar (and main) factor 
of the energy gap suppression is the phase jump $ which 
leads to the formation of a "potential well" around the 
junction having a width of the order of £o and contain- 
ing localized excitations with an energy \e\ < A (see Sec. 
3). In contrast to the ballistic case, the Josephson trans- 
port in a diffusive junction is performed not only by the 
states in the potential well, but by excitations within the 
whole energy region near the gap edge where the density 
of states differs significantly from the unperturbed value. 

II. EQUATIONS FOR GREEN'S FUNCTION OF A 
LOW-TRANSPARENT JOSEPHSON JUNCTION 

In order to calculate the density of states and equilib- 
rium supercurrent 

j = -v F v F D \ def (e)Tra z (g R Vg R -g A \7g A )(e) (7) 

we must solve equations for the matrix retarded (ad- 
vanced) Green's functions g R ' A (r,e) averaged over the 
ensemble of scatterers: 
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[a z e + Aexp(ia z x)i(Jy, g] =iDV(g\/g), (8) 

Here A and x are the modulus and phase of the order 
parameter and /o(e) = (l/2)(l + tanh(e/2T)) is the equi- 
librium distribution function. 

According to the normalization condition g 2 = 1 for 
the Green's function, the matrix g can be presented as 
g = cru, where er is the vector formed by Pauli matrices. 
Using the well-known relations (tra)(crb) — ab+icr[axb], 
[a z , er] = 2i [er x s] , where s is the unit vector of "isotopic 
spin" directed along the z-axis in the space of Pauli ma- 
trices, we can obtain from Eqs. (3) and (8) the follow- 
ing equations and the boundary conditions for the vector 
Green's function u: 

e[s xu}+ iA[x xu] = (L>/2)V[it x Vu], u 2 = 1, (9) 

€o[uxT7u]l,r = 2W[u l xu r ], (10) 

where \ = (sin;\;,cosX)0) is the symbolic vector of the 
order parameter phase. 

Singling out the component of the vector u along the 
direction s: u = su + iv (vs = 0), we project Eq. (9) 
onto the (a;, y)-plane in the space of Pauli matrices: 

ev - Aux = (iD/2)V(uVv - vVu), u 2 - v 2 = 1 (11) 

and introduce the unit vector ip = (simp, cos'0, 0) di- 
rected along v. v = tpv, where ip(r, e) is the phase of 
"anomalous" Green's function v (Vi/> = [tp x s]Vip). The 
obtained system of scalar equations is a possible repre- 
sentation of Usadell equations: 

ev-Aucos(ip-x) = ^D [V(uVv-vVu)-uv(Vtp) 2 ], (12) 

Aw sin(^ -x) = (j£>/2)V(v 2 VV>), (13) 

u 2 -v 2 = 1, (14) 
and its solutions determine the supercurrcnt 

/ + OC 
defo Im (v R ) 2 Vip R . (15) 
-oo 

Choosing the coordinate axis x orthogonally to the 
contact plane x = (x(+0) = -%(— 0) = $/2) and tak- 
ing into account the continuity of Green's function and 
antisymmetry of their derivatives, we can easily obtain 
from Eq. (10) the boundary conditions to Eqs. (12), (13) 
for x — > +0: 

£ {uVv - i>Vu)(0) = 4Wu(0)v(0) sin 2 tp(0), (16) 
C V^(0) = 2tUsin2V>(0). (17) 



Far away from the junction, the behavior of the order 
parameter and Green's function phases is described by 
linear asymptotic form corresponding to the given value 
of current 

x(+oo)=V>(+oo)=Xoc + 2 Ps x, Ps = (W/£,o)sm<S>, (18) 

i.e., of the superfluid momentum p s whose magnitude is 
determined in the main approximation by the condition 
of equality of the current Eq. (4) through the junction to 
its value j = 7re^F-Dp;5Atanh(A/2T) in the bulk of the 
metal. The Green's functions tend to their asymptotic 
values satisfying Eqs. (12)-(14) for ip — x an d Vw = 
Vu = 0. 

Using the parametrization u — coshf?,-u = sinh^, 
which takes into account the normalization condition Eq. 
(14), we can put in correspondence to the vector Green's 
function u the following geometrical image fl6|| . The unit 
vector u in a normal metal is directed along the isospin 
axis z (which corresponds to a purely electron or hole 
state of excitation of a Fermi gas), while in a supercon- 
ductor this vector is deflected from the axis through an 
imaginary angle iQ and turned around it through the az- 
imuthal angle tp. In the spatially homogeneous case, this 
angle obviously coincides with the phase of the order pa- 
rameter (tp = x)i an d the scalar Green's functions u and 
v are described by the formulas 

u R > A = cosh6 s = , £ u fl ^=sinh0 s , (19) 

\/(e± iO) 2 - A 2 ^ ' 

where ±i0 defines the position of singularities of the re- 
tarded (advanced) Green's function in the complex plane 
e, and the square root in Eq. (20) is defined so that 
u R ' A — > ±1 for e — * +oo. 

Eqs. (12)-(14) for Green's functions should be supple- 
mented by the self-consistency conditions for the modu- 
lus and phase of the order parameter: 

/+oo 
def a Re v R , (20) 
-OO 

/ + CO 
de/ Re« i? sin(^- X )=0, (21) 
-oo 

where A is the constant of superconducting interaction. 
Taking into account the current conservation law, Eqs. 
(13) and (21), it is convenient to calculate the value of 
current at the barrier (x — > +0) by expressing V"0(O) 
in Eq. (15) with the help of Eq. (17) through the phase 
jump 2-0(0): 

e f + °° 

j($) = --v F v F T / def Im (u fl (0)) 2 sin 2^ fl (0), (22) 

which allows us to single out explicitly the small parame- 
ter of the theory, i.e., the barrier transmissivity T. It can 
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easily be verified that in the main approximation using 
the unperturbed values of Green's function of Eq. (19) 
and phase -0(0) « x(0) = $/2, Eq. (22) leads to the 
result of Eq. (4). 

A simplifying factor in the case of a low transmissivity 
of the barrier is that the quantities ip — x an d Vi/> pro- 
portional to the current through the junction are small 
(see Eqs. (17) and (13)), and hence we can omit in Eq. 
(12) the terms quadratic in W and containing the phase 
gradients. Replacing ip(0) ~ x(0) = $/2 in the boundary 
conditions Eqs. (16), (17), to the same degree of accuracy, 
we obtain the equation and the boundary conditions for 
the parameter 9: 



esinhtf- A(x)cosh6» = (iD/2)V 2 



(23) 



£ V6>(0) = 2W sinh 26>(0) sin 2 $/2, 0(+oo) = 6 S . (24) 

Direct application of the perturbation theory to the so- 
lution of Eq. (23) (9{x) = 9 +Bi{x), A(x) = A + Ai(x)) 
leads to an expression for the correction 0\(x) containing 
nonintegrable singularities at the gap boundaries, and 
as a consequence, to the divergence of the corresponding 
correction to the Josephson current Eq. (4) . This is asso- 
ciated with the emergence of localized states of quasipar- 
ticles at a tunnel junction in the current-carrying state 
mentioned in Introduction and considered in the next 
section. 



III. LOCALIZED STATES AT A TUNNEL 
BARRIER 

It will be proved below that the depth of the "potential 
well" in the vicinity of the barrier is much larger than the 
scale of variation of the order parameter. Consequently, 
it is sufficient to confine an analysis of the behavior of the 
density of states to the model with a constant A, in which 
Eq. (23) has a simple solution describing the attenuation 
of perturbations of Green's functions at a distance ~ £o 
from the barrier: 



In the general case (for an arbitrary e), the solution of 
Eq. (27) can be obtained only numerically, but the pres- 
ence of the small parameter 7 in (26) and (27) makes it 
possible to apply the perturbation theory. Far away from 
the spectrum boundary, we can put 9(0) = 9 S on right- 
hand side of (26), which leads to the following expression 
for the correction to the density of states at the barrier: 



N(e,0) - N (e) = -2j Re 



isinh 6e sinh 20. 



(28) 



that becomes obviously inapplicable for |e| — > A, where 
\6 S \ — ► 00. In this region, we must apply the improved 
perturbation theory (IPT) by putting \z\, \z a \ ^> 1 for an 
arbitrary (not necessarily small) value of z — z s . This not 
only reduces the power of the general Eq. (27), but also 
allows us to write it in a universal form which does not 
contain the depairing parameter 7: 



(yVE-1) 2 



ly 



(29a) 
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= z/[3V2, £ = /3 2 (e-A)/A, /3=(2/ 7 ) 1 / 5 » 1, (29b) 



Relations Eq. (29) show that the increase in the den- 
sity of states is bounded by a quantity of the order of 
(3 ~ VF~ 2 / 5 as we approach the spectrum boundary. 
Thus, the range of applicability of the conventional per- 
turbation theory, Eq. (28), is determined by the condition 
(e — A)/ A 3> (3~ 2 and overlaps with the region of appli- 
cability (e - A)/A < 1 of the IPT. The boundary e* of 
the spectrum (the position of the bottom of the potential 
well), below which the density of states vanishes, corre- 
sponds to the emergence of a purely imaginary root of 
Eq. (29a) at the point E* = -(25/6)(2/3) 1 / 5 w -3.842: 



£*(<&) = A 



l-G [WsnS- 



C = 



25 



3-6V5 



i 5.824. (30) 



tanh 



6{x) 



tanh 



9(0) 



■ exp(— A; e |a;|), (25a) 



fc £ 2 = it; 2 sinh0 s , Re k e > 0. 



(25b) 



The dependence of the position of the spectrum bound- 
ary on the phase jump at the junction is illustrated by 
Fig. ^ in which a similar dependence of the position of the 
Andreev level Eq. (6) in a junction between pure super- 
conductors is shown for comparison. It should be noted 
that the scale of variation of e»(<i>) is much larger than 
the splitting of the Andreev level from the boundary of 
the continuous spectrum for the same barrier transmis- 
sivity. This is associated with the large value of the de- 
pairing parameter 7 in the diffusive junction as compared 
to the splitting parameter T of the Andreev level as well 
as with the large numerical value of the constant C defin- 
ing the shift of the spectrum boundary Eq. (30). Fig. |^ 
shows the results of numerical calculation of the density 
of states at the junction on the basis of the general for- 
mula Eq. (27) for different values of the depairing param- 
eter, which show that in addition of the root singularity 
2z 6 (z — z s ) 2 = ij 2 (z 2 — l)(z 4 — l) 2 , z s =exp0 s . (27) (~ \[i — e*) at the spectrum boundary, the quantity N(e) 



The quantity 9(0) satisfies the boundary condition fol- 
lowing from Eqs. (24) and (25): 

fc e £ sinh ° s = 7 sinh 20(0), 7 = ITsin 2 *<l, (26) 

which can be reduced to the eighth-power algebraic equa- 
tion in z = exp0(O): 
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has a "beak- type" root singularity for e — A. Its physical 
nature is associated with an infinite increase in the atten- 
uation length k^ 1 of the perturbation of Green's function 
in the bulk of the metal, Eq. (25), within the vicinity of 
the gap boundary. 



For e* < e < A, the density of states decreases expo- 
nentially with increasing distance from the junction (Fig. 
||), which corresponds qualitatively to the image of the 
potential well of depth A — and of width ~ £ with 
excitations localized in it. 




<J>/Jt 



FIG. 1. Phase dependence of the position of the bottom 
of the "potential well" e* , Eq. (30) , and the order parameter 
A(0), Eq. (48), in the vicinity of the tunnel junction (solid 
curves) at T = 0, W = 0.01 and fo/J = 5. The dashed curve 
shows for comparison the position of the Andreev level in 
a pure single-mode junction, Eq. (5), for the same barrier 
transmissivity. 
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FIG. 2. Dependence of the density of states iV(e, $, 0) at 
the tunnel junction on the energy of quasiparticles for various 
values of the depairing parameter 7 = IF sin 2 (<E>/2) (solid 
curves). The dashed curve shows the energy dependence of 
the unperturbed density of states No(e). 



FIG. 3. Energy dependence of the density of states 
N(e, 0) for 7 = 0.01 at different distances x from a diffusive 
tunnel junction: (curve 1), £0 (curve 2), and 5£o (curve 3). 

It is well known that the Josephson current is carried 
through a ballistic junction by localized excitations only 
and can be presented in the following form: 

»(*) 



j($) = -2e^^^tanh- 



2T 



(31) 



where the index n labels Andreev levels. At the same 
time, Eq. (22) for current expressed in the IPT approxi- 
mation in terms of the reduced variables of Eq. (29), 



*)« -1(A) tanh — sin § J — Im {y R f= j ($) 



(32) 



£,(*) 



shows that the charge transfer in a diffusive junction is 
performed not only by the states within the potential well 
(E < 0, e < A), but also by the excitations with energy 
e > A in the region e — A ~ A/3" 2 , where the den- 
sity of states differs significantly from the unperturbed 
value N (e). It should be noted in this connection that 
Argaman jlTj proposed an analog of Eq. (31) for a diffu- 
sive system, which can be obtained by the replacement 
of the energy e„($) of Andreev levels by the local value 
e(£, <E>, x) of the excitation energy for x = 0, which is adi- 
abatically deformed by supercurrent, using instead of the 
discrete number n the continuous variable 



de'N(e',<f>,x) 



(33) 
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viz., the number of states with an energy smaller than 
e (£ = 0(e 2 — A 2 )a/ e 2 — A 2 for a homogeneous super- 
conductor) pl| . One can assume that the contributions 
from the bound and delocalized states to the Josephson 
current are taken into account simultaneously by the for- 
mula 



-2evp 



(34) 



9$ 



2T 



which, however, leads to correct results only in the case of 
a homogeneous current-carrying state (where V% plays 
the role of $) or a wide S'iVS'-junction (with a width 
L £o of the normal layer) and is inapplicable for a 
narrow bridge and tunnel junction. Nevertheless, the 
consideration of the function e(£, x) is useful in these 
cases also since this allows us to visualize the variation 
of the energy distribution of quasiparticle states in the 
vicinity of the junction (Fig. ||). 
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FIG. 4. Lines corresponding to the number of states of 
quasiparticles £(e, a;) = const (Eq. (33)) for 7 = 0.01 in 
the vicinity of the junction. The dashed line shows the posi- 
tion of the bottom of the "potential well" (£ = 0, e = e*(<&)). 



IV. CURRENT-PHASE DEPENDENCE FOR A 
JUNCTION IN THE SECOND ORDER IN W 

Although the modified perturbation theory for Green's 
function in the energy representation described in the 
preceding section is the most physically obvious method 
operating with actual excitation energies, it leads to con- 
siderable formal difficulties in the calculation of correc- 
tions to the Josephson current Eq. (4). Indeed, it was 
shown in the previous section that the expression for j ($) 
calculated on the basis of the IPT for Green's functions, 
Eq. (32), coincides with Eq. (4) since the small IPT pa- 
rameter j3~ 2 cancels out as we go over to the reduced 



variables of Eq. (29). Thus, in order to calculate the cor- 
rections to Eq. (4) we are interested in, we must leave 
the approximation of Eq. (29) that describes the behav- 
ior of Green's functions correctly only in a narrow vicin- 
ity of singularity in the density of states. For this pur- 
pose, it is convenient to use the formalism of temperature 
Green's functions by going over from integration over en- 
ergy in Eqs. (20)-(22) to summation over the Matsubara 
frequencies u„ = 7rT(2n + 1), n = 0, ±1, ±2, . . .: 

j($) = -neu F v F rT Re w2 (°) sin 2^(0), (35) 

uj„>0 



A(x) = -2ttAT Y Im v(x) 



(36) 



UJ n >0 



and making the substitution e — > ico n in Eq. (23). This 
allows us to avoid divergences of the type of Eq. (28) in 
the perturbation theory which, unlike the IPT, makes it 
possible to take into account the coordinate dependence 
A(x). 

It is expedient to use as the main approximation in the 
asymptotic expansion = # + 0\ + • ■ ■ the "adiabatic" 
value of Green's function corresponding to the local value 
of A{x) = A + Ai(a;), (A x (oo) = 0): 



uq (x) = cosh 9q (x) 



u n (x) 



, vq (x) = sinh 6*o (x) ■ 



A 



uj n (x) 



(37) 



where Cb n (x) — y/ui n 2 + A 2 (a;). In this case, the correc- 
tion 0i (x) satisfies the nonhomogeneous equation 



V 2 1 - k 2 Jx = V 2 o , kl = 2Q n /D 



(38) 



with the boundary conditions V0i(+O) = 2IFsinh2# s x 
sin 2 < I > /2, 6*1(00) = 0, where cosh(9 s = u! n /u) n is the value 
of the Green's function far away from the junction with 
the unperturbed value of A, and Cj n — V^n 2 + A 2 . 

The self-consistency condition for Ai(x) following from 
Eq. (20), 



A i(l) T E 4" = - T E ^ Im g) (39) 



completes the system of equations for determining the 
corrections 0\ and Ai, whose solution in the Fourier rep- 
resentation has the form 



Ai(g) = ~8IFA-^Uin 2 $/2, 



(40a) 



lOJ r . 



e 1 (iu n ,q)=8WA- 



A(0) 



W„ q 2 + kl Z Q A(q) 



sin 2 $/2, (40b) 



A 2 



A{q) = A(0) + q 2 B(q), A(0) = 2vrT ]T — , (41a) 



UJ n >0 



6 



= ^ T E a 



w„>0 



+00 



o;3 q 2 + k 2 ' 



(41b) 



(^(iwn, x), Axfx)) =/ ^e^(^^„, g), Ai(?)) . 
J-00 27T 

As regards the correction to the asymptotic value Eq. 
(18) of the phase tp(x) of the Green's function, it is equal 
to zero in this approximation. In order to prove this, we 
introduce the quantity ip = ip — x <C 1, which, according 
to Eq. (13), obeys the equation 



VV - klip 



-V 2 



Xu 



(42) 



where xi — x( x ) ~ x(°°) is a correction to Eq. (18) local- 
ized near the junction. Taking into account the bound- 
ary condition Vyj(O) = — Vxi(O) following from Eqs. (17) 
and (18), we find that this equation has the simple so- 
lution <p(iuj n , q) — — g 2 Xi(g)/(g 2 + k 2 ) which leads, after 
the substitution into the self-consistency condition Eq. 
(21), to the homogeneous integral equation for xi(<?) : 



E 

UJ n >0 



A 



/•+00 



-00 



q 2 cosqx 

d( i 2TT2 Xi(g = °- ( 43 ) 
q- + kf 



The only nonsingular solution of Eq. (43) is Xi(<?) = 0> 
which proves the absence of a correction to the Joseph- 
son current due to the deviation of the behavior of the 
phases of the order parameter and Green's functions from 
the linear law Eq. (18). This result can be explained as 
follows. The correction Xi( x ) is obviously of the order 
of the small correction p s i(x) to the constant value p s 
of Eq. (18) in the vicinity of the junction, that ensures 
the conservation of the current upon a change in N(e) 
and A. Since the value of p s ~ W, the correction to this 
quantity, and hence Xi( x ) an d V have a higher order of 
smallness (~ W 2 ) than the corrections of the order of W 
we are interested in. 

Substituting Eqs. (40), (41) into Eq. (22), we obtain 
the required correction to the Josephson current: 



4T 



<Jj=j($)-jo($) = --r-/(A)sin$ £ Re 



w„>0 



A^ 



= -I(A)W Z(T) [ sin $ - - sin 2$ 



-f-OG 

Z(T) = ^^AA TJ24 Ij 



dk 



k 2 B(k) 



(44) 



, (45) 



where k^ = uj n /A, A(k) and B(k) are defined by Eqs. 
(41) upon the substitution fc w — > k u , and Wo and Ao are 
the values of W and A at T = 0. 



At low temperatures (T <C A), the summation over ui n 
in Eqs. (41) and (45) can be replaced by integration with 
respect to the continuous variable u>: 



A(0) = 1, B{k) 



tanh v dv 
k 2 + cosh v 




^ - 2\/l — k 2 arctan ^ 



1 -k 2 



k 2 



which leads to the following asymptotic value of the func- 
tion Z(T) for T -> 0: 



Z(T): 



8 



IT 2 , 



dk 



TTk 2 



. 2B 2 (k) 

(I + P) 9 / 4 l + k 2 B(k) 



<2.178. (46) 



In the vicinity of critical temperature (A <C T), the 
quantity A(0) sa 7C(3) A 2 /4?r 2 T 2 is small, and the main 
contribution to integral of Eq. (45) comes from the region 
of small wave vectors k ~ A/T corresponding to damping 
of perturbations at large distances of the order of £(T) cx 
(T c — T) -1 / 2 . This allows us to replace the function B(k) 
by its value ttA/4T for k = 0: 



B(0) dk 



(2n+l) 2 7 A(0) + k 2 B(0) 



2n, 



1 ttAq 
7C(3)T C 



5.099. 



(47) 



The results of numerical calculations of the Z(T) depen- 
dence within the entire temperature range < T < T c 
are presented in Fig. ||. 



N 




0.0 0.2 0.4 0.6 0.8 1.0 

T/T c 

FIG. 5. The function Z(T), Eq. (46), defining the temper- 
ature dependence of the ratio Sj($)/I(A), Eq. (44). 



7 



Similarly, we can calculate by using Eqs. (40) and (41) 
the asymptotic values of the correction Ai(0) to the un- 
perturbed value of the order parameter at the junction: 

^-=-a(T)Wosin 2 ^, a(0) = 3.037, a(T c ) = 5.782. (48) 
A 2 

The dependence of the order parameter A(0) on the 
phase jump at the junction at T = presented in Fig. 
1 shows that the main contribution to the energy gap 
suppression comes from the depairing mechanism con- 
sidered in Sec. 3, and the change in the order parameter 
is smaller than the variation of £*('&). 

The structure of the phase and temperature depen- 
dences of the correction to the Josephson current of Eq. 
(44) in a diffusive superconductor virtually coincide with 
expression Eq. (1) for a junction between pure metals 
except the following circumstance noted in Introduction: 
the parameter of the expansion of j($) in the transmis- 
sivity of the junction for I <C Co is not the tunneling 
probability T, but a considerably larger parameter W, 
Eq. (5). This allows one to observe higher harmonics of 
the current -phase dependence in diffusive tunnel junc- 
tion with a comparatively high resistance. Koops et al. 
plf apparently reported on the first experimental results 
in this field. 

The theory discussed above describes the current- 
phase dependence for a diffusive Josephson junction in 
the whole temperature range < T < T c except a nar- 
row neighborhood of T c in which A/T c ~ W (A/T c ~ T 
in a pure superconductor), and the magnitude of correc- 
tions Eqs. (44) and (1) becomes comparable with jo(&), 
while the correction Eq. (48) to A becomes of the order of 
its unperturbed value. This means that in the definition 
Eq. (5) of the parameter W near T c the coherence length 
Co (T) describing the characteristic scale of spatial varia- 
tions of Green's function and density of states should be 
replaced by the characteristic length £(T) of variation of 
the order parameter (healing length) in the Ginzburg- 
Landau theory, whose order of magnitude is the same 
as Co far away from T c . Taking into account the results 
of calculations of for a pure superconductor in the 
vicinity of T c §, we can obtain the following interpola- 
tion estimate of the effective transmissivity W suitable 
for any temperatures and mean free paths: 



(49) 



As we approach T c , the value of W increases infinitely, 
which is accompanied with a decrease in the phase jump 
for a given external current bounded by its critical value. 
Thus, in the ID geometry for an arbitrarily large normal 
resistance of the junction, there exists a narrow region 
near T c in which the phase difference of the order param- 
eter at the junction is small up to values of current of the 
order of the bulk critical current. 
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